Dependence of kinetic friction on velocity: Master equation approach 
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We investigate the velocity dependence of kinetic friction with a model which makes minimal 
assumptions on the actual mechanism of friction so that it can be applied at many scales provided 
the system involves multi-contact friction. Using a recently developed master equation approach 
we investigate the influence of two concurrent processes. First, at a nonzero temperature thermal 
fluctuations allow an activated breaking of contacts which are still below the threshold. As a result, 
the friction force monotonically increases with velocity. Second, the aging of contacts leads to a 
decrease of the friction force with velocity. Aging effects include two aspects: the delay in contact 
formation and aging of a contact itself, i.e., the change of its characteristics with the duration of 
stationary contact. All these processes are considered simultaneously with the master equation 
approach, giving a complete dependence of the kinetic friction force on the driving velocity and 
system temperature, provided the interface parameters are known. 

PACS numbers: 81.40.Pq; 46.55.+d; 61.72.Hh 



I. INTRODUCTION 

Almost three centuries ago Charles Coulomb (1736- 
1806) discovered that kinetic friction docs not depend on 
the sliding velocity [l|. Later, more careful experiments 
showed that this law is only approximately valid 043- 
Friction does depend on the sliding velocity, but this de- 
pendence is far from universal: some measurements find 
an increase when velocity increases, while others find a 
decay llOj] or even a more complex non-monotonous 
behavior A logarithmic dependence, often quoted, 
has been found for two extreme scales, friction at the tip 
of an AFM (see for instance M, fnHl3j ) or at the scale 
of a fault in the earth crust [l0|, but it is often only ap- 
proximate and observed in a fairly narrow velocity range. 
Therefore understanding the velocity dependence of ki- 
netic friction is still an open problem, and what makes it 
difficult is that several phenomena contribute, the ther- 
mal depining of contacts, their aging, and the delay in 
contact formation. 

In this study we investigate the velocity dependence 
of friction with a model that includes these three contri- 
butions and makes minimal assumptions on the actual 
mechanism of friction so that it can be applied at many 
scales provided the system involves multi-contact fric- 
tion. Our aim is to elucidate the respective role of these 
three contributions to the velocity dependence of friction 
and to provide analytical treatments in some limits, or 
simple numerical approaches that allow the investigation 
of a velocity range that may span many orders of magni- 
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tude. 

At the most fundamental level multi-contact friction 
can be described as resulting a succession of breaking and 
formation of local contacts which possess a distribution of 
breaking thresholds. This viewpoint was first applied to 
describe earthquakes [H|,[II| and then adopted to friction 
by Persson |16j |. 

We recently developed a master equation (ME) 
approach to describe the breaking and attachment 
events [13, EU • It splits the analysis in two independent 
parts: (i) the calculation of the friction force, given by the 
master equation provided the statistical properties of the 
contacts are known, and (ii) the study of the properties 
of the contact themselves, which is system dependent. 
This method is very general and allows us to calculate 
the velocity dependence of friction, which results from 
the interplay of two concurrent processes. First, at a 
nonzero temperature thermal fluctuations allow an ac- 
tivated breaking of contacts which are still below their 
mechanical breaking threshold. This phenomenon leads 
to a monotonic increase of the friction force F with the 
velocity v. Second, the aging of contacts [l9l - r2l1 | leads to 
a decrease of the friction force with velocity. It includes 
two processes: the delay in contact formation, i.e., time 
lag between contact breaking and re-making [lj|, |20|-[23| , 
and the aging of a contact itself, i.e., the change of its 
characteristics with the time of stationary contact. To 
incorporate the latter effect, the master equation must 
be completed by an equation for the evolution of static 
thresholds. 

In earlier studies [I?], EH we considered thermal and 
aging effects separately, to set up the method. How- 
ever, to relate the results to experiments, both con- 
tributions must be taken into account simultaneously. 
This is the aim of the present paper, which is orga- 
nized as follows. Section |TT] is a brief review of the mas- 
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ter equation approach. Section [TTT] discusses tempera- 
ture effects. Whereas our earlier work [l8[ focussed on 
time-dependent phenomena to analyze stick-slip, here we 
concentrate on the steady-state case (constant velocity). 
This allows us to proceed further and derive explicit ex- 
pressions for the influence of temperature. Then Sec. IIVI 
introduces the second effect, the aging of the contacts. 
It first summarizes the method introduced earlier and 
its main results [l8j], which only considered the T = 
case, and then studies the combined influence of aging 
and temperature fluctuations. Section fVl adds the influ- 
ence of the delay in contact formation after breaking, to 
get the full picture, allowing us to compute the velocity 
dependence of friction. Section [VTI discusses all those re- 
sults in the context of experimental data. The difficulty 
to apply the theory to actual experiments is to prop- 
erly assess the values of the parameters that enter in the 
theoretical expressions, and not simply try to fit experi- 
mental curves, which would not be very significant owing 
to the number of parameters which are involved. There- 
fore Sec. lVII focusses on this assessment. Finally, Sec. lVIII 
concludes the paper with a discussion of perspectives for 
its further development. 



II. MASTER EQUATION 

The earthquake (EQ) model is the most generic model 
for friction due to multiple contacts at an interface. The 
sliding interface is treated as a set of N c "contacts" which 
deform elastically with the average rigidity k. The con- 
tacts represent, for example, asperities for the interface of 
rough surfaces |24|. or patches of lubricant or its domains 
("solid islands" [25j |) in the case of lubricated friction. 
The ith contact connects the slider and the bottom sub- 
strate through a spring of elastic constant fcj. When the 
slider moves, the position of the contact point changes, 
and the contact's spring elongates or shortens, so that the 
slider experiences the force — F = fi from the inter- 
face, where fi — kiXi and Xi(t) is the spring length. The 
contacts are coupled frictionally to the slider. Namely, as 
long as the force \fi\ is below a certain threshold f S i (cor- 
responding to the onset of plastic flow of the entangled 
asperity, or to local shear-induced melting of the bound- 
ary lubrication layer), this contact i moves together with 
the slider. When the force exceeds the threshold, the con- 
tact breaks, and then re-attaches again in the unstressed 
state after some delay time r. Thus with every contact 
we may associate the threshold value f S i, which takes 
random values from a distribution P c {f) having a mean 
value f s . The spring constants are related to the thresh- 
old forces by the relationship fcj = fc (/si/is) 1 , because 
the value of the static threshold is proportional to the 
area Ai of the given contact, while the transverse rigid- 
ity fcj is proportional to contact's size, k{ oc \AAj' When 
a contact is formed again (re-attached to the slider), new 
values for its parameters have to be assigned. 

Rather than studying the evolution of the EQ model 



by numerical simulation it is possible to describe it ana- 
lytically [13, EH- Let P c {x) be the normalized probability 
distribution of values of the thresholds x S i at which con- 
tacts break; it is coupled with the distribution of thresh- 
old forces by the relationship P c (x) dx = P c (f) df. To de- 
scribe the evolution of the model, we introduce the distri- 
bution Q(x; X) of the stretchings Xi when the bottom of 
the solid block is at a position X. Let us consider a small 
displacement AX > of the bottom of the sliding block. 
It induces a variation of the stretching x.i of the contacts 
which has the same value AX for all contacts (here we 
neglect the elastic deformation of the block). The dis- 
placement X leads to three kinds of changes in the dis- 
tribution Q(x; X): first, there is a shift due to the global 
increase of the stretching of the asperities; second, some 
contacts break because their stretching exceeds the max- 
imum value that they can withstand; and third, those 
broken contacts form again, at a lower stretching, after 
a slip at the scale of the asperities, which locally reduces 
the tension within the corresponding asperities. These 
three contributions can be written as a master equation 
for Q(x;X): 



Q(x;X) = R(x)T(X), (1) 



where P(x) AX describes the fraction of contacts that 
break when the slider position changes from X to X + 
AX. At zero temperature P(x) is coupled with the 
threshold distribution P c {x) by the relationship fl7l. [iljj 

/>00 

P(x) = P c (x)/J c (x) , J c (x) = / d£P e (0 . (2) 

J X 

The function T(X) in Eq. ((T|) describes the contacts that 
form again after breaking, 



T(X) 



d£P(Z)Q(S;X) 



(3) 



(the delay time is neglected at this stage), and R(x) is 
the (normalized) distribution of stretchings for newborn 
contacts. Then, the friction force is given by 



/oo 
dxxQ(x;X). 
-OO 



(4) 



The evolution of the system in the quasi-static limit 
where inertia effects are neglected shows that, in the long 
term, the initial distribution approaches a stationary dis- 
tribution Q s (x) and the total force F becomes indepen- 
dent of X. This statement is valid for any distribution 
P c (x) except for the singular case of P c (x) = S(x — x s ). 

In the present work we concentrate on the steady state 
(smooth sliding). In what follows we use R(x) — 5(x) for 
simplicity. The steady-state solution of Eq. (TTJ) is 



Q(x)=Q(x)E P (x)/C[P) 



(5) 



where Q(x) is the Heaviside step function (Q(x) = 1 for 
x > and otherwise), Ep(x) = exp[— U(x)], U(x) — 
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So d t, p (0, and C[P] = J °° dxEp(x). Note also that, in 
the steady state, 



r = i/c[p] , 

because / °° d£P(£)E P (t) = J °° dU e~ u = 1. 



(6) 



The distribution Pr.(.f ) can be estimated for the contact 



of rough surfaces [18|, 124 
polycrystal substrates [18 
approximated by the function 



as well as for the contact of 
[26J: its general shape may be 



Pc(/)oc/ n exp(-///*), 



(7) 



where n > depends on the nature of the interface. 
Then, the distribution P c (x) can be related to the dis- 
tribution P c (fs) of the static friction force thresholds 
of the contacts. If a given contact has an area A, 
then it is characterized by the static friction threshold 
f s oc A and the (shear) elastic constant k cx VA (as- 
suming that the linear size of the contact and its height 
are of the same order of magnitude, see [l|| and Ap- 
pendix A in Ref. [HI)- The displacement threshold for 
the given contact is x s = f s /k, so that f s oc x s , or 
df s /dx s oc x s . Then, using P c (x s )dx s = P c (f s )df s , we 
obtain P c (x s ) oc x s P c [f s (x s )], or 



P c {x) oc x i+2n exp(-x'/xi) 



(8) 



where x* may be estimated from experiments as N c kx* ~ 
F a - In the SFA/B (surface force apparatus/balance) ex- 
periments, where the sliding surfaces are made of mica, 
the interface may be atomically flat over a macroscopic 
area. But even in this case the lubricant film cannot be 
ideally homogeneous throughout the whole contact area 
- it should be split into domains, e.g., with different 
orientation, because this will lower the system free en- 
ergy due to the increase of entropy. Domains of different 
orientations have different values for the thresholds / S j, 
i.e., they play the same role as asperities in the contact 
of rough surfaces. 

For the normalized distribution of static thresholds 
given by Eq. (|S} with n = 1, 



P c (x) = (2/x*)uV 



where u = x/x* 



(9) 



we can express the steady-state solution of the master 
equation analytically. In this case 

J c (x) = (1 + u 2 ) e~ u2 , (10) 
so that at zero temperature we have 

P(x) = (2/x*) U 3 /(i + « 2 ), (11) 



U(x) = ii 2 — ln(l +u 2 ) , 



E P (x) = J c (x) = (1 + u 2 )e 



(12) 



(13) 



C[P] = £*/Co, where Co = 



30? 



0.752, (14) 



Q{x) = (C /x*)(l + u 2 )e~ u , u>0 (15) 
and the kinetic friction is 

fk = F k / (N c k) = f k0 = C x* . (16) 

The ME formalism described above can be extended 
to take into account various generalizations of the EQ 
model, such as temperature effects and contact aging, 
which arc examined in the following sections. 



III. NONZERO TEMPERATURE 

Temperature effects enter in the ME formalism 
through their effect on the fraction of contacts that break 
per unit displacement of the sliding block, P(x), because 
thermal fluctuations allow an activated breaking of any 
contact which is still below the threshold [l6|, I181T20I-425 
For a sliding at velocity v so that X = vt, the thermally 
activated jumps can be incorporated in the master equa- 
tion, if we use, instead of the zero-temperature breaking 
fraction density P{x), an expression Pt{x) defined by 
(see [TJ) 



P T (x) = P{x) + H{x) 



(17) 



H(x) = - e kx2 ' 2kBT 



where the temperature contribution is given by 

d£P c (0e- ke/2kBT (18) 

v J x 
for "soft" contacts or by 

H(x) = £ fdtP c (0 (l - A\-M*(i-i)%/2k B T 

(19) 

in the case of "stiff" contacts which have a deep pinning 
potential so that their breaking only occurs with a sig- 
nificant probability when their stretching is close to the 
threshold. Here uj is the attempt frequency of contact 
breaking, u ~ 10 s accor ding to Refs. QSHlj]. 

For concreteness, in what follows we assume that the 
contacts are soft, Eq. (p~8|) . and we select n = 1 in Eq. 
so that P c (x) is given by Eq. ©. 

At a nonzero temperature the total rate of contact 
breaking, Eq. (TTT)) . is equal to Pt (x) = P(x)+(w/v) h(x), 
where for the soft contacts 



[) (1 + 6) 2 



with 



b(T) 



kxj 
2k B T 



(20) 



(21) 
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The condition 6=1 defines a crossover temperature 

k B T* = -kx\. (22) 

Then, a straightforward integration gives the function 
U T (x) = /* d£iV(0 = t/"(iB) + A[/(s), where 



At/ (a:) =5 («,T) erf («) - S x (T) u e~ u ' 



with 



and 



6/3) 



5o(u ' T) = ^7TT6F 



_ Co (1 + 6) 

5l(T) -T(TT673) 



(23) 



(24) 



(25) 



The coefhcient Si(T) weakly changes with temperature 
from 5i(0) = 2/0F ~ 1.128 to Si(oo) = C /2 w 0.376. 
On the other hand, the coefficient So(v,T) determines 
whether the effect of temperature is essential or not. The 
temperature-induced breaking of contacts is essential at 
low driving velocities only, when Sq{v, T) 3> 1. Thus, the 
equation So(v*,T) = 1 defines the crossover velocity: 



w.(T) 



ux* [1 + 6(T)/3] 
Co [1 + 6(T)P 



(26) 



We see that monotonically increases with temperature 
as v*(T) « 0.443 wx»T/T* at T < T* and approaches the 
maximal value w 1.33 ukc* at T > T». 

Then, E PT {x) = e ~ UT ^ = (1 + u 2 ) e^e"^*), and 
we can find the kinetic friction force: 



fk(v,T) 



I (IxxEpt(x) I I dxEpT 



Or). (27) 



At a low driving velocity, v -C w * , we may put At/ (or) 
Sgli, where 



v{l + b) 2 



and Eq. fl27]) leads to 

f k ~X,/S 2 = {v/L0){l + bf 



(28) 



(29) 



A linear dependence of the kinetic friction on the driving 
velocity at low velocities corresponds to the creep mo- 
tion due to temperature activated breaking of contacts 
and was predicted in several earlier studies [ljl [2l| - |23l | , 
although our approach allowed us to derive it rigorously. 
The dependence Q29p may be interpreted as an effective 
"viscosity" of the confined interface: 



fk 1 , , 
r\ = — = - 1 

V U! 



kxi y 

2k B T J 



(30) 




FIG. 1: (Color online): The friction force fk as a function of 
the driving velocity. Dash-dotted blue line shows the low-t) 
approximation (|29[) . dashed red line shows the high-u approx- 
imation (I3ip . and dotted magenta line shows a logarithmic 
fitting, k = 1, CjJ = 1, x* = 1, feT = 1. 



At a high velocity, v 3> v», when e Ac/ ( ;r ) «l-5; 



^(^u) 2 , we obtain C 
that 



2 u- 



3v^/4 - S* 2 + (5V?/16)5| so 



fk ~ fkO (l — C1S2 + C2S2) , 



(31) 



where Ci = 50F/8-C o « 0.356 and C 2 = 16/9vr-l/3« 
0.233. Equation (f3"TT) agrees qualitatively with that found 
by Persson [16( in the case of 6 ^> 1. 

Approximate expressions (|29[ l3Tj) together with the 
numerical integration of Eq. ()27|) are presented in Fig. [T] 
Also we showed a logarithmic fitting which operates in 
a narrow interval of velocities near the crossover ve- 
locity only. Persson showed that the logarithmic 
dependence may be obtained analytically, only if the 
P c {x) distribution has a sharp cutoff at some x = x s 
as, e.g., in simplified versions of the EQ model with 
P c (x) = 5(x - x s ). 

Although we cannot obtain analytical results for the 
stiff contacts, we calculated the fk(v) dependences nu- 
merically (see Fig. [5]), which shows that the effect remains 
qualitatively the same. 



IV. AGING OF CONTACTS 



The aging of contacts was considered in our work 18 1 
where, however, we ignored the temperature-induced 
breaking of contacts discussed above in Sec. IIIII When 
aging is taken into account the master equation for 
Q(x, X) must be completed by an equation for the evolu- 
tion of P c (x), which in turns affects P(x). Let the new- 
born contacts be characterized by a distribution P C i(x), 
while at t — ¥ 00, due to aging the distribution P c (x) ap- 
proaches a final distribution P c f{x). If we assume that 
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0.01 




FIG. 2: (Color online): The friction force ft as a function of 
the driving velocity v for soft (dotted) and stiff (solid curves) 
contacts at low temperature ksT = 0.01 (blue) and high tem- 
perature ksT = 1 (black), fe = 1, w = 1, as* = 1. 



the evolution of P c (x) corresponds to a stochastic pro- 
cess, then it should be described by a Smoluchowsky 
equation 

dP - - 8 ( 3 \ 

-af = DLP °> where is + W 

in which the "diffusion" parameter D describes the rate 
of aging, B{x) — dU(x)/dx, and the "potential" U(x) 
determines the final distribution, P c j{x) oc exp [ — f7 (a;)] , 
so that we can write 



B(x) = - [dP cf (x)/dx] /P cf {x) 



(33) 



However, because the contacts continuously break and 
form again when the substrate moves, this introduces two 
extra contributions in the equation determining 8P c /dX 
in addition to the pure aging effect described by Eq. ([32]) : 
a term P(x; X)Q(x; X) takes into account the contacts 
that break, while their reappearance with the threshold 
distribution P C i (x) gives rise to the second extra term in 
the equation. Thus, the evolution of P c is described by 
the equation 



dP c (x;X) 
dX 



-D v LP c {x; X) + P(x; X) Q(x; X) 
= P ci (x)r(X), (34) 



where D v — D/v, and v = dX(t)/dt is the driving veloc- 
ity. Finally, we come to the set of equations (jT] [3l [34]) . 
For the steady-state regime, Eq. (jM)) reduces to 

D V C[P] LP c (x) = P(x) E P {x) - P ci {x) , (35) 

where we used Eqs. © and ©. Taking also into account 
the identity P(x) Ep{x) = P c (x) [la ], we finally come to 
the equation 



FIG. 3: (Color online): The kinetic friction force fk as a 
function of the driving velocity v for different values of the 
aging rate: D = oo (red dotted), 1 (blue short-dashed), 10 _1 
(black solid), 10 -2 (magenta dash-dotted), 10 -3 (blue short- 
dotted), and D = (red dotted curve), k = 1, w = 1, fcsT = 
1; the initial and final P c (x) distributions are given by Eq. ([9]) 
with x,i = 0.1 and x t f = 1 correspondingly. Dashed curves 
shows the dependences at T — 0. 



It was shown [l8| that the kinetic friction monoton- 
ically decreases with the driving velocity as Fk(v) — 
-Ffc(O) oc —v/D in the low- velocity limit and Fk(v) — 
Fk(oo) oc D/v in the high- velocity case. One may ex- 
pect that at low velocities this decreasing will compen- 
sate the friction increasing due to temperature induced 
jumps. The problem, however, is more involved. 

When the temperature effects are incorporated, 
Eq. (|36| for the function P c (x) in the steady state takes 
the form 



D v C[Fr] LP c {x) = P c (x) - P a {x) . 



(37) 



Numerical solutions of Eq. (13"T|) are presented in Figs. [3] 
and 3^,: the initial increase of the kinetic friction F with 
the driving velocity v due to the temperature activated 
breaking of contacts is followed by the decrease of F due 
to contacts aging. FigureHb shows also the dependence of 
the effective "viscosity" rj = fk/v on the driving velocity 
v. It is constant at low velocity and then decreases; the 
latter may be approximately fitted by a power law r](v) oc 
v~ a with the exponent a changing from 1.5 to 1 as D 
decreases. 

Using the definition (|3"12|) of the operator L and Eq. (|3"3"1) 
for the function B(x), the l.h.s. of Eq. (I3T1) may be rewrit- 
ten as 

while the r.h.s. of Eq. (|37p may be presented as 



D V C[P] LP c (x) = P c {x) - P cl {x) . 



(36) 



P c (x)-P ci {x) 



dx 



[J c {x) - Jd] , 



(39) 



6 



o.i 



0.01 

1 



H 



0.1 



0.01 




D = 1 
-£>= 10 
■D= 10 
•D=10 




7 = 0.31 v 
-- ^ = 0.105 v" 
,) 0.08 v 



0.1 



FIG. 4: (Color online): (a) The same as in Fig. [3] in log- 
log scale, (b) The dependences of the effective "viscosity" 
r\ — fk/v on the velocity v (dashed lines show power-law fits). 



where J cl {x) = f™d£P C i(Q. Using Eqs. ([38) and (|39|) . 
we can find the first integral of Eq. (|37j) : 



FIG. 5: (Color online): The kinetic friction force fk as a 
function of the driving velocity v for different values of the 
aging rate D = 1 (magenta dashed), 10 _1 (blue solid), 10 -2 
(black dash-dotted) and 10 -3 (red dotted curve) as compared 
with approximate expressions (thin solid curves), k = f, cj = 
1, ksT = 1; the initial and final P c (x) distributions are given 
by Eq. ([9]) with x ti = O.f and x,f = 1 correspondingly. 



where Ui — x/x„ and 

1 1 



D v C[P T }P cf (x)±(^- 
dx \P c f{x) 



9 2 2 

■^if x */ 



(44) 



Jdix) — J c (x) . (40) Then, taking the corresponding integrals, we obtain to 



first order in 



Integration of Eq. (|40]) leads to an integral equation for 
the function P c (x): 



C ~ x^i (Cq 1 — S*2 + C12S3 



P c (x)=P cf (x)\l + 



DC[P T ] J P cf (0 

(41) 

Substituting J c (x) « J c f(x) = f^°d{;P c f(t;) into the 
r.h.s. of Eq. f4"Tj) , one may analytically find the low- 
velocity behavior of the kinetic friction, for example, the 
decrease of fk with v for T = 0. At a nonzero tem- 
perature, however, aging does not affect the low-velocity 
behavior $Z§§ and only reduces the interval of velocities 
where Eq. (J2U) is valid, as demonstrated in Figs. |3] and 2] 
Indeed, at v — > and T > the main contribution to 
Pt{x) comes from the function H(x) oc lj/v in Eq. (|17p. 
which only weakly depends on P c {x). 

The limit D v = D/v — > may be studied with the 
help of Eq. (|3"Tf by substituting P c (x) « P C i{x) into its 
left-hand side. For the function © , this approach gives 



and 



where 



~C 



2 r 



1 — (Co + C\)S2 + C12S2 



S 3 (v) = 4D v C[P]/x 2 if . 



(45) 



(46) 



(47) 



Pc(x) ~ Pd(x) 

and 



if 



"2 ) 



(42) 



(43) 



C12 ~ 1.324 and C12 ~ 1.789 are numerical constants. A 
comparison of the exact and approximate expressions is 
shown in Fig. [5] 



V. DELAY IN CONTACT FORMATION 

Finally, let us take into account the delay in contact 
formation following the work of Schallamach 22] . Let r 
be the delay time, N be the total number of contacts, N c 
be the number of coupled (pinned) contacts, and Nf — 
N—N c be the number of detached (sliding) contacts. The 
fraction of contacts that detach per unit displacement of 
the sliding block is T(v,T) = J dx P(x) Q(x), i.e., when 
the slider shifts by A A, the number of detached con- 
tacts changes by A c rAA, so that Nf — TvtN c . Us- 
ing N c + N f = N, we obtain N c = N/(l + Tvt) and 
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FIG. 6: (Color online): The kinetic friction force fk as a 
function of the driving velocity v for different values of the 
delay time: Vd — oo (black solid), 10 (red dot-dashed), 1 
(magenta dotted) and 0.1 (blue dashed curve) for 73 = 0.1 
and fcflT = 1 (other parameters as in Fig.[3|. Inset shows the 
same in log-log scale. 



N f = NTvt/(1 + Tvt). 
Vd = x/t, we can write 



If we define x = 1/T and 



N 



1 + v/vd 



and N 



f 



1 + v/vd 



(48) 



The coupled contacts produce the force fk defined 
above by the steady-state solution of the master equa- 
tion. The combined dependence which incorporates tem- 
perature effects, aging and delay in contact formation, is 
shown in Fig. [5] for different values of the parameter Vd- 

However, above we assumed that the sliding contacts 
experience zero friction, while these contacts may expe- 
rience a viscous friction force fi = rjiv, where rji corre- 
sponds to the (bulk) viscosity of the liquid lubricant. In 
this case the kinetic friction should be additionally mul- 
tiplied by a factor fj(v) = 1 + v /v v Vd, where v v — fk/r]i 
(Vf) 3> Vd). Such a correction may be expected at huge 
velocities only, e.g., for v ~ 1 m/s. In this case the 
function Fk(v), after decreasing, reaches a minimum at a 
velocity v ps (v^Vd) 1 ^ 2 , and then increases according to 
a law fk(v) oc rjiv. Note that the viscous friction which 
comes from the excitation of phonons in the substrates, 
as shown in MD simulation j7j, may also depend on the 
velocity, e.g., as r\i cx v 4 . 



VI. MAKING THE LINK WITH EXPERIMENTS 

For a real system, the results presented in the previ- 
ous sections allow the calculation of the kinetic friction 
force Fk(v,T) provided the parameters of the model are 
known. In this section section we examine how they can 
be evaluated from experiments. 

The contact parameters k and w may be estimated 
with the help of elastic theory [28| . Let us assume that 



a contact has a cylinder shape of height h (the thickness 
of the interface) and radius r c , so that it is character- 
ized by the section A; = irr 2 , the (geometrical) inertial 
momentum / = irrf /A, a mass density p and a Young 
modulus E. If the cylinder foot is fixed and a force A/ 
is applied to its top, the latter will be shifted on the dis- 
tance Ax = Afh 3 /3EI (the problem of bending pivot, 
see Sec. 20, example 3 in Ref. [28j). Thus, the effec- 
tive elastic constant of the contact is k = Af/Ax = 
3EI/h 3 . The minimal frequency of bending vibration of 
the pivot with one fixed end and one free end, is given 
by w « (3.52/ h 2 ){EI / pA,) 1 / 2 (see Sec. 25, example 6 in 
Ref. (Hi). 

Next, let a be the average distance between the con- 
tacts, so that the total area of the interface is A = Na 2 , 
and introduce the dimensionless parameter 7 C = r c /a 
(7 C < 0.5). The threshold distance x* may be estimated 
as follows. At the beginning, when all contacts are in 
the unstressed state, the maximal force the slider may 
sustain is equal to s» Nkx* (this force corresponds to 
the first large stick spike in the F(t) dependence at the 
beginning of stick-slip motion at low driving). Thus, we 
obtain that kx* ~ a 2 <r*, where <r* = F^/A is the maximal 
shear stress. 

Let us consider a contact of two rough surfaces and 
assume that a = h = r r . Then we obtain 



L76y/g7j 

u s» 

for the attempt frequency, 

k = (37r/4)£r c 
for the contact elasticity, and 

x * — t c o * 1 1 k 



(49) 



(50) 



(51) 



for the threshold distance. For steel substrates we may 
take p = 10 4 kg/m 3 for the mass density, E = 2 x 
10 11 N/m 2 for the Young modulus, and a c = 10 9 N/m 2 
for the plasticity threshold. Assuming that er* = a c 
and r c ps 1 pm, we hnd that to pa 7.9 x 10 9 s _1 , 
k pa 4.7 x 10 5 N/m, pa 2.1 x 10~ 9 m, b pa 2.7 x 10 8 
for room temperature (i.e., b 3> 1), so that the crossover 
velocity is quite low, v* sa 0.03 pm/s. 

Now let us consider a lubricated system, e.g., the one 
with a few OMCTS layers as studied by Klein [29[ and 
Bureau [27|, and assume that the lubricant consists of 
solidified islands which melt under stress as proposed by 
Persson [25|. In this case, instead of using the Young 
modulus, let us assume that x* — r c ; this allows us to 
find the parameter EI = ah 3 a*/3j c . Then, the elastic 
constant is k — aa*/j c , the attempt frequency is 



u pa 1.15 \J a*/ahp-f 3 , 
the parameter b is given by 

b sa 7 c a 3 a*/2fc B T , 



(52) 



(53) 



8 



CU in 4 




dyldt (s~ ) 

FIG. 7: The shear stress a as a function of the shear rate 7 
for different values of the OMCTS film thickness from m = 2 
to 6 monolayers. Inset: the crossover velocity v* as a function 
of the number of layers. 



and in the case of b ^> 1 the crossover velocity is 

pa j c u;a/3Cob ps fc B T '/ 'y/a^hpay^ . (54) 

For a four-layer OMCTS film [29J one may take p = 
956 kg/m 3 , h ps 3.5 x 10~ 9 m, F* ps 2 x lfT 5 N and A ps 
10~ 10 m 2 so that cr» « 2x 10 5 Pa. Assuming 7 C = 0.5 and 
a pa 1 /im, we obtain for room temperature, fcsT = 4 x 
1(T 21 J, that w sa 8 x 10 8 s" 1 and b pa 1.25 x 10 7 , i.e. this 
system is in the low-temperature limit too, although the 
crossover velocity is much higher than for rough surfaces, 
w* pa 16 jtmi/s. 

Moreover, we may calculate the dependence fk(v) for 
different thicknesses of the lubricant film. If the film con- 
sists of n\ layers, then the film thickness is h = n\d, where 
d pa 8.75 A is the diameter of the OMCTS molecule. Let 
us assume that the maximal shear stress exponentially 
decreases with the number of layers according to the re- 
sults of MD simulation [7|, a* = coe"^™ 1 , where ~ 1 
is a numerical constant. Taking <tq = 4 x 10 6 N/m 2 and 
/3 = 1.5, we obtain the dependences of the shear stress 
u = Fk/A on the shear rate 7 = v/h shown in Fig. 
which may be compared with the experimental depen- 
dences (Fig. 2a) of Bureau [27 [ . 

Note that our approach may overestimate the value of 
the crossover velocity v*. First, the crossover will occur 
earlier if the delay and/or aging effects play a signifi- 
cant role. Besides, at low temperatures the stiff contacts 
lead to higher "viscosity" and lower values of u* than 
the soft contacts considered above (see Fig. [5]). Second, 
we completely ignored the elastic interaction between the 
contacts. If the latter would be incorporated, a break- 
ing of one contact may stimulate neighboring contacts to 
break as well, i.e., the value of the parameter a should 
describe such a cooperative "contact" size which may be 
much larger than those of individual ones. 



Giving a quantitative evaluation of the influence of ag- 
ing on the velocity dependence of the friction coefficient 
is harder than for the temperature dependence due to 
insufficient experimental data. Aging appears to cause a 
decrease of friction as velocity increases, and thus, when 
such a behavior is observed experimentally 0, [l3[ it can 
be considered as a strong indication of the presence of 
aging. Our analysis indicates that the combined effect 
of temperature and aging leads to a maximum in the 
friction coefficient versus velocity. Therefore, when ag- 
ing is manifested by a decreasing friction versus velocity, 
extending the experiments to lower velocities and tem- 
peratures might detect the maximum and thus provide 
some quantitative data to evaluate the aging parameters. 

Although the aim of our work was to find the depen- 
dence of the kinetic friction on the driving velocity, our 
approach allows us to find the dependence on tempera- 
ture as well. However, the behavior of a real tribologi- 
cal system is more involved, because all parameters may 
depend on temperature T in a general case. For exam- 
ple, the delay time r may exponentially depend on T 
if the formation of a new contact is an activated pro- 
cess 22[; the same may be true for the aging rate D. 



In this case one may obtain a nonmonotonic tempera- 
ture dependence of friction with, e.g., a peak at cryogenic 
temperatures [2lj . 



VII. CONCLUSION 

In this study we determined the dependence of the ki- 
netic friction force in the smooth sliding regime on the 
driving velocity. In a general case, the friction linearly 
increases with the velocity (this creep motion may be in- 
terpreted as an effective "viscosity" of the confined film) , 
passes through a maximum and then decreases due to 
delay/aging effects. The decay may be followed by a new 
growth in friction in the case of liquid lubricant. Esti- 
mation showed that for the contact of rough surfaces, 
the initial growth of friction should occur at quite low 
velocities, v <C 0.1 /im/s, so that for typical velocities 
the friction is independent on velocity in agreement with 
the Coulomb law. However, for the case of lubricated 
friction with a thin lubricant film which solidifies due 
to compression, the fk(v) dependence is essential, and 
the linear dependence may stay valid up to velocities 
v ~ 10 4- 10 3 /im/s. At higher velocities the growth satu- 
rates and the fk{v) dependence may be fitted by a loga- 
rithmic law. The latter velocity interval is narrow if the 
distribution of static thresholds is wide; the logarithmic 
law may be found analytically for a wide interval of ve- 
locities when the thresholds are approximately identical, 
i.e., for the singular distribution P c (x) — S(x — x s ). 

We emphasize that our approach is only valid for a 
system with many contacts, for example, N > 20 at 
least [2lJ. When the contact is due to a single atom 
as it may occur in the AFM/FFM devices, the friction 
can be accurately described by the Prandtl-Tomlinson 



9 



model and should follow the logarithmic fk(v) depen- 
dence, f k (v) cx (lnw/w ) 2 / 3 :',() . But if the AFM/FFM 
tip is not too sharp so that the contact is due to more 
than one atom, the logarithmic dependence is only ap- 
proximate and, moreover, for some systems the friction 
may decrease with the velocity which has to be attributed 
to the aging /delay effects d,S,|3l|. 

In this work we had in mind that contacts correspond 
to real asperities in the case of the contact of rough sur- 
faces or to "solid islands" for the lubricated interface. 
However, the ME approach also operates when the con- 
tact is due to long molecules which are attached by their 
ends to both substrates. Such a system was first studied 
by Schallamach [22j and then further investigated by Fil- 
ippov et al. [13], Srinivasan and Walcott (23[, and Barel 
et al. [21]. Note that when all molecules are identical, 
they are characterized by the same static threshold, i.e., 
this system is close to the singular one, where the loga- 
rithmic fk(v) dependence has to have a wide interval of 
operation. 

Finally, let us discuss restrictions of our approach. 
First of all, we assumed the somehow idealized case of 
wearless friction; wearing may mask the predicted de- 



pendences. Besides, the interface is heated during slid- 
ing; this effect is hard to describe analytically as well 
as to control experimentally. Then, we did not estimated 
the delay /aging parameters; moreover, these parameters, 
e.g., the delay time r, may depend on the driving velocity 
v. Besides, we assumed the simplest mechanism of aging 
described by the Smoluchowsky equation, while the real 
situation may be more involved, e.g. , it may correspond 
to the Lifshitz-Slozov mechanism 18]. Also, we assumed 
that the reformed contacts appear in the unstressed state, 
R(x) = S(x) in Eq. (JXJ) , which may not be the case in real 
systems. 

The most important issue, however, is the incorpora- 
tion of the elastic interaction between the contacts as 
well as elastic deformation of substrates at sliding. This 
point certainly deserves a detailed investigation and is 
the topic of our future work. 
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